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SYMBOLS 


van Driest damping coefficient 

modeling constant in term representing turbulence dissipation, 
equations (2) and (14) 

specific heat at constant pressure 

turbulence function, equation (8) 

specific kinetic energy of turbulence in in.compressible flow 
function of turbulence Reynolds number, equation (10) 

■v 

mass-averaged stagnation enthalpy of fluid 
mass-averaged fluctuating static enthalpy 
Karman constant in mixing length formulation 

modeling constaiit in texmi representing the redistribution between 
Reynolds stresses 

mass- averaged specific turbulence kinetic energy, pu^UjJ/p 
length scale of turbulence 
mixing length 

turbulence function, equation (7) 
mean local static pressure 
fluctuating local static pressure 

instantaneous heat flux in ith direction by molecular means 

Reynolds number of turbulence 

empirical modeling constant, equation (13) 

correlation length used by Glushko 

instantaneous rate of strain tensor 

mass-averaged mean rate of strain tensor 

mass-averaged fluctuating rate of strain tensor 

instantaneous temperature 



T mass- averaged mean temperature 

T" fluctuating temperature in mass-averaged coordinates 

u, Uj^ mean velocity in streamwise direction 
u^ mass-averaged mean velocity in ith direction 

u' fluctuating velocity in streamwise direction 

u^ . fluctuating velocity in ith direction 
u^ mass-averaged fluctuating velocity in ith direction 

V, U2 mean velocity normal to a surface 
V* fluctuating velocity normal to a surface 


X, y I 

xi, X2I 
a 

r 

5. 

«ik 

e 

< 


X 

V 


p 


< >k 
(-) 


coordinates in directions along and normal to a surface 

empirical modeling constant, equation (12) 

modeling constant related to a turbulent Prandtl number 

viscous layer thickness 

Kronecker delta, i = k, = Ij i ^ k, = 0 

ratio of eddy diffusivity to kinematic viscosity 

modeling constant in term representing turbulence dissipation, 
equation (2) 

empirical modeling constant,, equations (8) and (15) 

✓ 

kinematic viscosity on incompressible fluid 
eddy diffusivity 

mean density in a compressible fluid 
density of an incompressible fluid 

instantaneous stress caused by nwlecular motions (viscosity) 
modeling constant in compressibility term, equation (66) 
partial differentiation in kth direction 
time-averaged quantity 
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A ONE-EQUATION MODEL OF TURBULENCE FOR USE WITH 


THE COMPRESSIBLE NAVIER-STOKES EQUATIONS 
Morris W. Rubes In 
Ames Research Center 


SUMMARY 


The Glushko one-equation model of turbulence is extended to compressible 
flows without boundary-layer approximations.) 

INTRODUCTION 


Recent advances in both computer technology and numerical analysis permit 
the evaltjation of compressible flow fields by numerically solving the complete 
Navier-Stokes equations, including the energy equation for evaluating variable 
fluid properties (e.g., ref. 1). Thus, flow fields containing separated 
regions, which hitherto had to be treated by iterative procedures that coupled 
essentially viscous- and inviscid-flow regions, can be treated as a single 
field and solved with reasonable computer economy (refs. 2-5). A difficulty 
arises, however, when these techniques are applied to turbulent flows through 
the Reynolds stress formalism. For separated flows caused by adverse pressure 
gradients, the comparativeiy large pressure gradients, both along and normal 
to body surfaces, and the presence of the separated- flow region in close 
proximity to a surface are factors outside the bounds of previous turbulence 
modeling experience. Nonetheless, mixing length models, expressed by empirical 
formulas evaluated from experiments on attached turbulent boundary layers in 
equilibrium in the Clauser sense (ref. 6) , have been tried in the numerical 
solutions and the results have been compared with experimental data on turbu- 
lent boundary layers separated by either incident shock waves or deflected 
surfaces such as ramps (refs. 3 and 4). Although these computations yield 
results for the flow-field quantities — surface skin friction and heat trans- 
fer that have the general character of the data — their quantitative differ- 
ences are rather large and in certain regions there are significant qualitative 
differences as well. To improve the computations, intuitive relaxation models 
were tried that use empirical formulas for mixing length and/or eddy viscosity 
lODdlfled by a rate equation based on the degree of departure from local equi- 
librium. Although the predictions of the pressure distributions were improved 
for some experimental cases, the resulting skin- friction values generally were 
not significantly better than the original models (ref. 7) . An early version 
of a two-equation model of turbulence was also tried (ref. 3) and the results 
were not any better than those resulting from the empirical models. 

The failure' of the early two-equation turbulence model to improve the 
calculations over the results from the empirical algebraic models illustrates 
that more generality and complexity in turbulence msdeling do not necessarily 



assure greater accuracy in a particular application. It was decided therefore 
to try a one-equation model for use with the compressible Navier-Stokes equa- 
tions, but one that contains more of the physics of turbulence than the relaxa- 
tion models and that is not as complex as a two-equation model. The incompres- 
sible fluid model of Glushko (ref, 8)' was adopted for this purpose because of 
its rather direct formulation, its lack of arbitrary damping functions, and 
its rather good representation of flat-plate boundary-layer flows. Also, this 
model has been carefully compared with a well-documented incompressible 
boundary- layer flow in an adverse pressure gradient (ref. 9); modeling modifi- 
cations based on the data are indicated there. The same authors also show 
(ref. 10), that data obtained using the Glushko model agree fairly well with 
Other boundairy-layer data obtained under other arbitrary pressure distributions . 

The present report suggests the modeling modifications necessary to 
extend the Glushko model to compressible flow within the full Navier-Stokes 
equations. Implied in this extension is the assumption of the adequacy of the 
empirical length scale inherent in the model to represent separated flows near 
the surface. This report first presents an outline of the Glushko model for 
incompressible boundary-layer flow. Then the compressible analog of the model 
is derived in the absence of the boundary-layer approximation. Compressibility 
introduces new terms modeled in two different ways to show the arbitrariness 
inherent in the modeling process. The tentative model, subject to future test 
against experiment, is then summarized. 


GLUSHKO MODEL 


The Glushko turbulence model for an incompressible, turbulent boundary- 
layer flow is used to guide the development of a one-*equation model of turbu- 
lence applicable to the Navier-Stokes equations for co mpressible flows. 

Glushko derived a model for the Reynolds stress, u'v', in the x-direction 
momentum equation for an incompressible turbulent boundary- layer flow. The 
Reynolds stress was shown to be expressible by the product of an eddy viscosity 
and the strain of the mean flow, 8u/3y, in what is usually called the logarith- 
mic region of the boundary layer. This eddy viscosity simplification was not 
assumed a priori , but resulted directly from the Reynolds stress equations in 
an analysis (ref. 11) using the boundary-layer assumptions and neglecting the 
convection and diffusion of the Reynolds stresses. The principal turbulence 
mechanisms used were the production, dissipation, and redistribution of the 
individual components of the Reynolds stresses. For the later mechanisms, 
Glushko used Rotta's models (ref. 12), namely, 

p'<'^2,k + '^k.i^ “ -^1 T (e-^ ■ I ^ik^®) 

for the redistribution of turbulence between the components of the Reynolds 
stress and 
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for the dissipation. 

In addition to establishing an eddy viscosity form for the Reynolds shear 
stress, Glushko's analysis also demonstrated that the eddy viscosity was 
expressible in terns of a local turbulence Reynolds number defined as 


Glushko ide ntified the length scale . L, with half the distance, S, where 
the correlation UjCxj, X2)uj(x2, Xo + S) vanishes. Although this definition 
of L is quite arbitrary in the modeling equations, L always occurs in a 
product with an empirical constant so that there is no loss of generality. 
Finally, Glushko did not use the relationship between e and r found 
directly from the analysis of the Reynolds stress equations , but used it to 
justify the use of an empirically developed relationship between these quanti- 
ties that presumably corrects for the previous neglect of the convection and 
diffusion of the individual components of the Reynolds stresses. 

The Glusko model for a two-dimensional incompressible, turbulent boundary 
layer is 


Continuity of mass: 


Suj 9 u 
_ ■' 0 
3 X 2 3 X 2 


Momentum balance in X2 direction: 
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£ 3X2 3X2 


[•Sr] 


Kinetic energy of turbulence: 


— 3 e , _ 3 e , •. 

“1 3^7 "2 ^ “ 


ST 


where 


M = 1 + e(r) 


D = 1 + e(Xr) 
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Note that the argument of e, the ratio of eddy diffusivity to kinematic vis- 
cosity, in the definition of D has been altered by the empirical constant \ 
to account for differences in the transport of momentum and kinetic energy. 

This can be considered as equivalent to introducing a Prandtl number for the 
transport of turbulence kinetic energy. The ratio of eddy diffusivity to kine- 
matic viscosity is expressed as 


e = H(r)ar 

where H(r) is expressed empirically as 


(9) 



'r/Vo . 

0 5 < 0,75 \ 

H(r) = - 

r/r^ - (r/r^ - 0.75)2 ^ 

0.75 < r/r^ < 1,25 



1.25 < r/r^ < » J 


( 10 ) 


The length scale is also given by an empirical expression based 
boundary- layer correlation length data, namely, 

on flat-plate. 


, 0 < X 2/6 < 0.23 


L/5 = 

< (x2/d + 0.37)/2.61 , 0.23 ^ x^/d 4 0.57 

1(1.48 - X2/d)/2.52 , 0.57 < X2/d < 1.48, 

( 11 ) 

With the empirical 

constants. 

. 


a » 0.2 

( 12 ) 


rp = 110 

(13) 


C = 3.93 

(14) 


X = 0.4 

(15) 


Equations (3) through (11) form a closed system that can be used to calculate 
incompressible turbulent boundary-layer flow. It is emphasized again that e, 
the specific turbulence kinetic energy, is found from a differential equation 
involving the production, dissipation, and diffusion of turbulent kinetic 
energy so that the turbulence need not be in equilibrium with the meari flow. 
This dynamic characteristic is introduced back into the mean-flow equations 
through the Reynolds number, r. 


COMPRESSIBLE ANALOG TO GLUSHKO MODEL 


To extend the foregoing model of turbulence to the compressible Navier- 
Stokes eqiaations, the use of mass-averaged dependent variables (ref. 13) is 
adopted. The mean field equations, written in abbreviated differential form, 
are 
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Continuity of mass: 


p,^. + (pUj),j = 0 


(16) 


Conservation of momentum in ith direction: 


(pUj^),j, + [pUjUj^ + ° 


(17) 


Conservation of thermal energy: 
(pH - P) 


' *f pu. 
,t 2 


puVuj 


H + + pu^'h" - Ui (Tij - pu^u|.') - It ^ 


= 0 


(18) 


'J>2 


The last term in equation (18) is usually neglected since it is smaller than 
the thermal dissipation term to it s left . The terms that must be evaluated 
then to close these equations are pujh" and pu^Uj. 


Consistent with the Glushko model, an eddy viscosity concept is adopted. 
For nonboundary- layer flows, however, an additional assumption that the eddy 
viscosity is a scalar quantity is necessary. With these assumptions, the com- 
ponents of the Reynolds stress tensor are written: 


■pUj^Uj - 


^€(ui^j + Uj^i - I 6ijUk,k) - I 6ijpk 


(19) 


in the notation of reference 13. The heat flux vector can be expressed as 

^ ^ fi.j (20) 

or, in keeping with the way Glushko handled the transport of turbulence 
kinetic energy. 


puJh" » -iie(rf)hj (21) 

where T is a new universal constant related to the turbulent Prandtl number, 
Pr^ and f is defined later in equation (35). 

From equation (19), the effective total shear stress, laminar plus turbu- 
lent, can be written as 

fij - puVuj = p(l + £)S^j - -| ^ijPk 

where the mean strain of the flow is 

®ij “ '^i,j '^j,i “ 1 ^ij'^k,k 


( 22 ) 


(23) 
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With equations (19) through (23), eqxiations (16) through (18) are closed once 
e is defined. 

To define e (the ratio of eddy diffusivity to kinematic viscosity) when 
the fluid is compressible, one must do the following; 

(1) Show that past experience makes it reasonable to use equation (11) 
for the length scale in compressible flow. 

(2) Show that s is a function of a turbulence Reynolds number in com- 
pressible flow and show what form this Reynolds number takes. 

(3) Model the turbulence kinetic energy equation to account for effects 
of compressibility. 


Length Scales in Compressible Flow 


For equilibrium turbulent boundary layers, it has been found that com- 
pressibility has little direct effect on the length scales of turbulence. For 
example, the early mixing length theories for a boundary layer on a flat plate 
(e.g., refs. 14 and 15) and the more recent computer codes for calculating the 
boundary-layer equations on two-dimensional bodies with mild pressure gradients 
(ref. 16) have shown good agreement with experimental data up to Mach numbers 
of 6 while treating the mixing length as being essentially independent of 
density. This independence can be illustrated quite clearly through a con- 
sideration of the turbulence model used in reference 16 where the eddy diffus- 
ivity is expressed in terms of an inner and outer part of the boundary layer. 

In either part of the boundary layer, reference 16 uses 


V 


t 



In the inner part, the mixing length scale is 


^inner “ " ®^P<“y/A)] 


(24) 


(25) 


with 


26 )J_ 

A = — (26) 

on a flat plate with zero surface mass transfer and 

K = 0.4 . (27) 

the usual Karman constant value. In reference 16, local properties are used 
as the reference fluid properties in the definition of A in equation (27) . 
This causes the value of A to depend somewhat on the density or compressi- 
bility of the fluid within the sublayer adjacent to the surface. Outside the 


6 



sublayer, however, y > A so that the exponential term in equation (25) becomes 
small and effects of compressibility tend to vanish from dinner* Thus, over 
most of the inner region of the boundary layer, the length scale is independent 
of the fluid density. 

In the outer part of the boundary layer, reference 16 uses the Clauser 
form of eddy diffuslvity: 



(28) 


Note that the density does not appear explicitly in this expression. Since the 
eddy diffusivity in the outer part of the boxmdary layer is independent of 
density, then, together with equation (24) , . the suggestion is that the length 
scale in the outer part of the boundary layer is also independent of the den- 
sity of the fluid, other than through the effect of compressibility on the 
boundary-layer thickness 6. This is seen more directly in the model used in 
reference 17, where outside the sublayer the mixing length is expressed as 



(29) 


so as to be independent of density. 


From the above discussion of mixing length used in programs that accu- 
rately predict near-equilibrium turbulent boundary-layer behavior, it can be 
concluded that the leJngth scale L in equation (11) will apply to compressible 
flows as well. Its application to separated flows is less clear and modifica- 
tions may be needed under those conditions. . 


Turbulence Reynolds Number in Gompressible Flow 

In a compressible turbulent boundary layer, the total shear stress 
expressed in equation (22) reduces to 

T = Tj2 - pu"u^ = iad + e) = pv(l + s) (30) • 

The mixing length expression, equation (24) , can be written in terms of the 
ratio of eddy diffusivity to the kinematic viscosity as 

ly 

Equations (30) and (31) combine to yield 

e(l + s) = (32) 

v2 
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In the fully turbulent portion of the boundary layer, e » 1 and 


_T 

P 


pu"v" 

p 


( 33 ) 


so that under these conditions 


/-(pu"v'Vp)& 

V 


(34) 


Thus, e identifies with a turbulent Reynolds number based on local fluid 
properties and the local correlation of the fluctuating velocities. This sug- 
gests, by analogy, that the Glushko' bkj del can be extended to compressible 
flows by use of equation (9), but with the turbulence Reynolds number f 
redefined as 


r 



(35) 


where ic is the mass-averaged kinetic energy of turbulence and v is the 
local kinematic viscosity. 


Second-Order Modeling of Turbulence in Kinetic Energy Equation 

The local mass-averaged specific kinetic energy of turbulence is repre- 
sented in reference 13 by 


(pk)^t + “ -<P'^i^j)"l,j- HP^i + P“i,i+K'^ij\j^ 

(36) 


The first term on the right In equation (36) can be related through equa- 
tion (19) to the mean velocity field and the compressible extension of equa- 
tion (9), namely, e = H(f)ar. The second term, having divergence form, repre- 
sents a diffusion of the instantaneous turbulence kinetic energy by the 
turbulence itself. The third and fourth terms .representing pressure and veloc- 
ity correlations and the fifth and sixth terms containing molecular-shear and 
velocity correlations must be modeled to close equation (36). 


Models of pressure and velocity correlations . - When the pressure in the 
abovementioned terms is represented as the sum of a mean value plus a fluc- 
tuating quantity, p = p + p' , the pressure and velocity correlation terms 
become 




(37) 


For inconq)ressible flow, both u^ = 0 and “ 0 so that only the second 

term on the right remains. This term, having a divergence form, is identified 
with the diffusion of turbulent pressure fluctuations by the turbulence and is 
usually grouped with the second term on the right in equation (36) . The other 
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terms in equation (37) require models that converge to zero when the flow is 
inGompressible. Thus, grouping these terms yields 


I’V + (uiP’),i * -u^p,i + (38) 

To evaluate the quantity u^ in the' first term on the right in equa- 
tion (38), it can be assumed first that the total temperature of the fluid is 
constant within a field of turbulent eddies: 


Tj. = T + * const 


(39) 


When T is replaced by a mean and fluctuating mass-averaged quantity (ref. 13), 


and 


Then 


T = f + T" 


“k ■ \ + >v 


1! 




2^p 


On averaging, equation (42) becomes 




2c 


Cp 2Cp 


Equation (43) subtracted from (42) yields 


ii„!t 


Cp • K 1C- ZCp 2Cp 


X" - T" + — (u," - u,") + = 

which is satisfied by 

% « 

0 


n tr'Vu," 

T" + — uj; + = 

Cp ^ 2Cp 


(40) 


(41) 


(42) 


(43) 


(44) 


(45) 


Since the mean velocity in at least one direction will be much larger than the 
fluctuating velocity in that direction, the above equation can be well repre- 
sented by its leading terms: 


T" = - ^ 


(46) 
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At this point, it is usual to argue that pressure fluctuations are of a 
smaller order than either density or temperature fluctuations , . and this assump- 
tion, together with the perfect gas equation, permits one to identify the tem- 
perature fluctuation with a density fluctuation in equation (46) , thus relating 
density and velocity fluctuations. Since this assumption is suspect at high 
Mach numbers, an alternative assumption will be proposed, namely, that the gas 
behaves in a poly tropic manner. Then 


= P* _ n T" 
P ~“p ~(n-l) X 


(47) 


where n is the polytropic coefficient (n = 0, isobaric; n = 1, isothermal; 
n= Cp/c^, isentropie, etc.)* 


From equations (46) and (47), the density fluctuation is related to the 
velocity fluctuation as 


P’ 


P ~ If 

(n - l)Cpf K 


,(48) 


With equation (48) , 


p’uV 


uV = - 


(n - Dcpf V^i 


(49) 


If the third-order correlations between the density and velocity are neglected, 

1 


uV = - 


(n - l)CpT 


- f ®ik^) 


(50) 


Since the value of n is not known, it, in effect, becomes a modeling coeffi- 
cient and equation (50) can be used to define the first term in equation (38). 
Note that n » for incompressible flow (see eq. (47)) and u^ ^ 0 as 
required. Under boundary- layer assumptions, the only appropriate value for 
the index is i = 1 and, consequently, equation (49) reduces to 


ti, = 


u^u”u” 

(n -^ l)cpf 


= const 


u^k 


(n - l)Cpf 


= const 


(Y - 1) “1*^ 
(n - 1) ;2 


(51) 


Since k/a^ is a measure of the square of the turbulence Mach number, u" is 
expected to be small if n takes on a value of about 1.2 as expected in non- 
isentropic flows. This latter form is identical with the Alber and Wilcox 
representation (ref. 18) , except that the n must be less than unity to have 
the same sign as in their work. 


An alternative representation of u^ that does not require the assump- 
tion of constant total temperature within the turbulence is 


u^ = 


p’uV 

X 


• -r. ^ T”u** 

(n - 1)T ^ '"i 


(52) 


10 



If density fluctuations are small compared to the mean density, then equa- 
tipn (52) can be rewritten as 


or 


uV = - 

X 


(n - Dc^Tp ^ (n 


Z , A) ijie(rf)h • 

Dp Y 


( 53 ) 


^ (Y - 1) ve(rf) 

^ - 1 ) 



(54) 


In the absence of an application of either equation (50) or (54) to predict 
several sets of experimental data, it is not clear at this point which is the 
better model, although equation (54) has the advantage of being simpler and not 
requiring the assumption of a constant total temperature within an eddy. This 
latter assumption turns out to be rather significant. For example, in boundary- 
lay er flow and with, say, n > 1, equation (51) would assign a positive value to 
u”* On the other hand, for the same conditions, equation (54) can assign 
either positive or negative values, depending on whether the flow was over a 
cooled or a heated surface. This illustrates how tenuous modeling assumptions 
are until established by comparison with data. 


With the assumption that state variables within an eddy behave in a poly- 
tropic fashion, the third term on the right in equation (38) is expressed as 


The instantaneous continuity equation 


(55) 


p + (pu.) . = 0 
> ^ J > J 


(56) 


when written in terms of mean and fluctuating quantities is 

p + p'^ + (pu. + puV + p'u, + p'u") ^ = 0 
3 J 3 3 >3 


(57) 


When equation (16) is subtracted from this equation, the instantaneous fluc- 
tuations of density and velocity are related by 


p’ + (puV + p'u, + p*uV) , = 0 (58) 

If the mean flow varies slowly over the scale of an eddy, and only the linear 
terms in the fluctuations are retained, equation (58) reduces to its highest 
order terms: 


p' + u.p' . + pu',' . = 0 
,t j ,j j,3 

With equation (59) multiplied by p* and then time-averaged. 


( 59 ) 
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The time- averaged quantities now scale with the dimensions of the mean flow 
field and can be compared with the mean continuity equation in similar form: 

p + up , + pu. . = 0 (61) 

> ^ J > J J 5 J 


If the rms density fluctuation intensity is called 6: 

<f) - 6 (62) 


then equation (60) becomes 


B^Cp + 

9 ^ 


u.p .) 

J .3 


+ p6(e 




+ u , 6 . ) + p'u'j 


3 .3 


3.3 


(63) 


It is now assumed that the intensity of turbulence changes relatively slowly 
along the streamlines in contrast to across the streamlines, so that the second 
term in equation (63) can be neglected. With equation (61), the resulting 
equation reduces to 

p^ . = &^pu. . (64) 

3.3 3.3 


The variation in turbulence intensity 6 from streamline to streamline, can 
be accounted for by setting 6^ proportional to the local kinetic energy of 
the turbulence. Usually, u^ > U 2 » even in separated flows, so that equa- 
tion (48) permits one to write 


, '"I k 

6“^ = const * ~ 

3 . 3 

which, together with equations (55) and (64), gives 

r, 2 . 

i k ... 


u 


i,i 


P'UV . = ?P — ^ U 


(65) 


( 66 ) 


where ^ is a modeling coefficient that includes the polytropic coefficient. 

Models of molecular-shear and velocity correlations .- The molecular- shear 
and velocity correlation terms in eqtiation (36) provide for the dissipation 
and diffusion, by molecular processes, of the specific turbulence kinetic 
energy. To demonstrate the assumptions underlying the modeling of some ele- 
ments of these terms, one must first expand the form of the terms to reveal 
their components. 



( 67 ) 


At an instant of time, the local molecular shear is expressed as 


T . . = pS . . 

11 IJ 

where y and are the instantaneous viscosity and strain. Specifically, 


’ij 


“l.j 3«ij"k,k 


( 68 ) 


When each variable is expressed as the sum of a mass-averaged mean quantity 
plus a fluctuation, equation (67) becomes 

Tij = yS^j + yS'^^ + y"S^^ + y"S”j (69) 

With equation (69), the molecular-shear and velocity correlation terms in 
equation (36) can be written as 


MSV . 


The terms containing y" are neglected'. The justification for this is that 
the third term should be negligibly small compared to the first because y" < y 
and the correlation coefficient between ul[ and y^j is much less than 
unity. By similar reasoning, the fourth term is small compared to the second. 
When the remaining terms are expanded, 


MSV = uVCyS..) . + (yuVS':,) . - yS'.'.u'.' . 
1 13 ,3 1 13 .3 13 1.3 




The first, fourth, and fifth terms then vanish in incompressible flow since 
u^ = 0 and u{J^k = 0. For compressible flow, the following arguments can be 
used to justify neglecting the fourth and fifth terms. Since they contain the 
second derivatives of fluctuating quantities before averaging, the smallest 
eddies should contribute most to these quantities. The smallest eddies, how- 
ever, tend to be isotropic (ref. 19), and, if the variation of p over the 
dimensions of these small eddies is neglect ed, s light extens ions of the 
methods in reference 20 (ch. 3)., show that uVu!! -• s = uVuV a a vanish for 
isotropxG turbulence xn a compressxble fluxd. 


The final expression then for the molecular-shear and velocity correla- 
tion term is 


MSV * uVCyS..) , + (yk ) 

1 13 .3 .3 .3 


- yuV ,uV , 
1.3 1.3 


(72) 


In keeping with the Glushko model, the third term representing the dissipation 
of the specific kinetic energy of the turbulence is modeled as 
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(73) 


iiuy ^ ^ = yC -5- [1 + e(Xr)] 

The in the first term can be taken from equation (54) , 

^fc>del of third— order correlation term .-^ The divergence terms on the right 
in equation (36) are grouped and inadeled as follows: 


(pu"k) + (u"pV) = Iu"(p' + puyu")] . 

J > J X ,X 3 X/ X. ,3 


« -ME 


(74) 


SUMMARY OF ONE-EQUATION TURBULENCE MODEL IN COMPRESSIBLE FLOW 


The system of equations that represents a one-equation model of compres- 
sible turbulent flow is as follows: 


Continuity of mass (see eq, (16)): 


p + (pu ) . = 0 

> ^ J > J 


(75) 


Conservation of momentum in ith direction (see eqs. (17), (22), and (23)) 


■ ° 


(76) 


where 


"^ij “ ‘^^i'^j “ e(r)]S_ - ^ 6^^pk 


(77) 


and 


S..=u. .+u. .--:r 6.,iL , 
ij i.J 3 xj^,k 

Conservation of thermal energy (see eqs. (18) and (21)): 


(pH - p) + [pu.H + q. + pu''h" - u. (t.. - puVuV)] = 0 

» t 3 J J X X3 X 3 , 3 


(78) 


(79) 


where 


q^ + PUj 


— —n I ^ 

j" --“[pr 


+ e(rr)Ki . 


(80) 


Specific turbulence kinetic-energy equation (see eqs. (36), (39), (54), (74), 
and (77)): 


(pk) + (pu k) . = -puyu'/u, , - (pu k) .-(uVp') . + [PV+ MSV+ (u'.'p ’ ) ] (81) 

> ^ J 9 J 4-3X53 3 > J J 9 j J 9 J 
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Cu"p' + pu''k) a -jie(Xr)k 
3 3 >3 >3 

PV + (“jP' ) J + p’u^ ^ 

ITT (y - 1) oecrf) 

"i ■ (n - 1 ) 52 

Ul^ k 

P u" = CP — ~2 u. ^ 
a‘^ J > j 


MSV = u"(qs ) . + (Pk ,) . - pc ;% [1 + e(Xr)] 
■*■ *3 »3 »3 »3 L 


These equations combine to yield 


(pk) ^ + (Pu^k) . = ye(r) S u. . - -r 5 pku + {1i[l + e(Ar)]k } 

>*- 3 *3 ^ 9 J ^ 9 J 9 3 93 

- $ “j.J + - W ^ [1 + (87) 

where n, C, and F are new modeling coefficients. Additional modeling quan- 
tities are 


e * H(r)ar 


r/fo , 


H(r) = <r/ro - (f/Pg - 0.75)' 


0 < v/Tq < 0.75 
0.75 < r/rg < 1.25 


1.25 ^ r/% < 


x,/5 , 


0 < x^/iS < 0.23 


L/6 “ < (X2/<S + 0.37)/2.61 , 0.23 < x^^/S < 0.57 

1(1.48 - x«/5)/2.52 , 0.57 < x«/5 < 1.48 
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where 


a = 0.2 I 

r^=110 I 

C = 3.93 
X = 0.4 

/ 

CONCLUDING REMARKS 


(92) 


The set of equations suimnarized in the preceding section is closed when 
the modeling parameters n, and T are established. To do this, one must 
eventually compare the predictions given by the model 'with the data from a 
series of experiments that covers a sizable range of Mach and Reynolds numbers. 
At this point, all that can be done is^ to offer suggestions for the first set 
of trial values of these parameters. For example, in nonisentropic inviscid 
flows, values of n — 1.2 seem to fit much of the data. Co nsistent with this 
value of n, Y = 1.4 and for the relationship be tween u”^ ^ that 

occurs in an equilibrium boundary layer, that is, u*’^ ^ 8/9 k, 5 ~ 8/11. The 
quantity T is identified with the turbulent Prandtl number, P = 1/Pr^, for 
which r = 1.1 is a reasonable first trial value. 
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